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Abstract 

We develop a model for gels and entangled solutions of semiflexible biopoly- 
mers such as F-actin. Such networks play a crucial structural role in the 
cytoskeleton of cells. We show that the rheologic properties of these networks 
can result from nonclassical rubber elasticity. This model can explain a num- 
ber of elastic properties of such networks in vitro, including the concentration 
dependence of the storage modulus and yield strain. 

PACS numbers: 61.25.H, 83.80.L, 87.10, 87.45 



1 



Typeset using REVTeX 



F.C. MacKintosh, J. Kas, and P.A. Janmey 



A variety of semiflexible biopolymers and protein filaments affect cell structure and 
function. The most prevalent of these in eucaryotic cells is actin, which forms the cytoskeletal 
rim [1,2]. This actin cortex is a polymer gel that provides mechanical stability to cells, and 
plays a key role in cell motion. Networks of actin and other protein filaments in vitro have 
been the subject of considerable recent interest [2-10], not only because of their structural role 
in cells, but also because of unusual viscoelastic properties of these networks. Such protein 
filaments as actin are novel in that they form viscoelastic networks, in which a <C £ ^ £ p , 
where a is the size of a monomer, £ is the characteristic "mesh" size of the network, and l v 
is the persistence length of a chain. In the case of actin, £ and £ p are of order 1/xm. This, for 
instance, has permitted direct visualization of polymer dynamics such as reptation [11,12] 
by optical microscopy [10]. Insight into the control of viscoelasticity in networks of both 
natural and synthetic semiflexible polymers in this intermediate regime is also important for 
the design of biocompatible materials. For instance, aqueous gels of stiff protein filaments 
or biocompatible polymers have both structural and pharmaceutical applications. However, 
neither models of flexible-chain solutions nor models of rigid-rod networks [11,12] are directly 
applicable to such systems. Here, we report a model for the elasticity of semiflexible polymer 
networks that can account for many of the observed properties of such networks in vitro. 

Solutions of actin filaments in vitro exhibit a polydisperse length distribution of about 
2 fj,m to 70 /im in length, with a mean length of 22 /xm [13]. Between 36 //g/ml and 2 
mg/ml F-actin forms entangled solutions. Many of the properties that are important for the 
function of the actin cortex appear to arise from essential differences of F-actin networks 
from gels and concentrated solutions of flexible polymer chains. We show that the rather 
large elastic moduli may not be generated by a classical entropic rubber elasticity, since the 
persistence length £ p of the filaments is comparable to or larger than the characteristic mesh 
size of the network £, as illustrated in Fig. |TJ This semiflexible nature of the filaments can 
also explain both qualitatively and quantitatively a number of elastic properties of actin gels 
and solutions. We shall focus primarily on actin networks, although our model is applicable 
to other semiflexible polymers at intermediate concentrations. 

Polymer gels and concentrated solutions are characterized by entangled chains. Depend- 
ing on preparation, these chains can be chemically crosslinked. For such a gel, the density of 
chemical crosslinks can be controlled by preparation. This density then determines the aver- 
age distance between points along the chain that are effectively fixed by the bonds with the 
surrounding elastic network of other chains. For an entangled solution, on the other hand, 
the viscoelastic properties depend on transient entanglements of an individual chain with 
its neighbors [11,12]. Despite the transient nature of these entanglements, over intermediate 
time scales of interest, the effect is much the same as that of chemical crosslinks, although 
the effective degree of entanglement or the average length L e between entanglements is more 
subtle [14]. This intermediate regime is the "rubber plateau", for which the solution behaves 
as an elastic solid. It is this regime that we address below. 

Although many viscoelastic properties of actin and other biopolymer networks resemble 
those of high molecular weight solutions of flexible polymer chains, the rubber plateau regime 
exhibits novel behavior. Actin solutions, for instance, exhibit relatively high plateau moduli, 
of order 100 Pa or higher for actin monomer concentrations of order 1 mg/ml (i.e., for volume 
fractions of order 0.1%) [9]. The plateau modulus of actin networks also exhibits significant 
strain hardening for modest strains. A rather small linear regime is observed: e.g., in many 



2 



F.C. MacKintosh, J. Kas, and P.A. Janmey 



cases they have a threshold strain as low as 10%, beyond which they lose their mechanical 
integrity. In the case of actin, this maximal strain also weakly decreases with increasing 
actin concentration [15]. As we show, this is a direct consequence of the intrinsic bending 
rigidity of biopolymers such as actin, and is direct evidence of the inapplicability of the 
freely-jointed chain model for the concentrations of interest [16-18]. 

We develop a model for densely crosslinked actin gels and entangled solutions, in which 
the elastic properties arise from chains that are very nearly straight between entanglements, 
as illustrated in Fig. [I]. As we shall focus on the elastic rubber plateau modulus, we shall 
not distinguish between crosslinked gels and entangled solutions, except in so far as the 
entanglement lengths may differ. We show that for an entangled solution, the plateau 
modulus scales with concentration ca as G' ~ c 1 ^ 5 . As shown in Fig. ^| this is consistent 
with the measurements to date of the concentration dependence of G' in the range of 0.3 — 
2.0 mg/ml [5]. For densely crosslinked gels, however, a somewhat stronger, G' ~ c^ 2 , 
dependence is predicted. 

In our model for the linear elasticity in the plateau regime, we consider an ensemble of 
chain segments of length L e (either between crosslinks or entanglement points), which are 
embedded in a continuous medium that undergoes a uniform shear deformation characterized 
by angle 9. The elastic response of the network results from the tension in such chain 
segments as a function of the extension, L — Lq, where Lq is the relaxed length. When a 
semiflexible chain segment is stretched by a tension r, the energy per unit length of the 
chain depends on two effects: the bending of the chain, and the work of contracting against 
the applied tension. The energy per unit length can be written [20] 

H= l -K(v*uf+ l -r(yu)\ (1) 

where k is the chain bending modulus, and u(x) describes the (transverse) deviation of the 
chain away from a straight conformation along the x axis, n is related to the persistence 
length of the chain i v (the length over which the chain appears straight in the presence of 
thermal undulations) by £ p ~ n/{kT). We let denote the full contour length of the 
chain (i.e., for k = oo or r = oo). We neglect the possibility of "internal" stretching of 
the chain: i.e., the chain is assumed to have no longitudinal compliance. Thus, for fixed 
contour length, — L ~ | / dx(Vu) 2 . At a given temperature and for a given tension 
r, the transverse thermal fluctuations of u determine the equilibrium length L. The chain 
conformation can be described by the Fourier series u{x) = J2 q u q sin(gx), where we include 
wavevectors q = j, . . . consistent with fixed ends of the chain segment. For the harmonic 

energy of Eq. ([]]), the mean square amplitudes (u 2 ^ can be calculated from the equipartition 
theorem, with the result that 



L^-L-kT^—^—, (2) 
nq 1 + r 

where we have included both transverse polarizations of u. To linear order in applied ten- 
sion r, the average end-to-end distance of the chain segment is L ~ — kTL 2 / (6k) + 
kT L A / (90 k 2 )t. The second term represents the equilibrium contraction of the end-to-end 
distance at finite temperature. The last term gives the linear relationship between the ap- 
plied tension and extension 5L of the chain segment beyond its relaxed length. For a small 
extension, the tension is given by [21] 
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T ~wr> 5L - < 3 » 

The above results for the behavior of individual chains can be used to estimate first the 
maximum shear strain # max that a network can withstand. This will, in general, decrease 
with increasing concentration, since the entanglement length will then decrease. This means 
that the fraction of the excess chain length in the form of thermal undulations decreases, 
and hence there is less chain available to "pull out" under the applied stress. More precisely, 
the relative extension of a segment of length L e between entanglements is proportional to 
the strain 9: 5L ~ 9L e . Considering the total excess length — Lq above, the maximum 
strain for chain segments of length L e is given by 8 max ~ kTL e /n. Thus, the maximum 
strain is predicted to depend linearly on L e . Furthermore, this maximal strain decreases 
with increasing chain stiffness (for the same entanglement length L e ). This is consistent 
with the observation that the yield strain does indeed increase with increasing flexibility of 
the network: networks of ADP actin, ATP actin, and vimentin show such a trend [17]. 

For the modulus, G', we use the relation above for the tension on an individual chain 
segment as a function of the shear strain in the linear regime. For a network, we consider 
such a chain segment of length L e that is deformed by an amount given by SL ~ 9L e . Within 
the linear regime (i.e., for small strain), the tension in the chain segment is thus given by Eq. 
(0). Solutions and gels are characterized by a mesh size £ that describes the average spacing 
between chains or the size of voids between filaments. Along a plane parallel to the shear, 
there are l/£ 2 chains per unit area [12]. The stress a is therefore given by a ~ k 2 / (kTC, 2 Ll)8 
in the linear regime. Thus the modulus scales as 

G' ~ ^C 2 Lf. (4) 

This is in contrast with the behavior of gels and networks of flexible chains, for which 
G' ~ kT/t; 3 [11]. 

Both the entanglement length L e and the mesh size £ decrease with increasing concen- 
tration of chains, although, unlike concentrated solutions of flexible chains, the scaling of 
these quantities with concentration need not be the the same when L e ?Z £ [22]. The char- 
acteristic mesh size £ for a network of stiff chains is given by £ ~ 1/ ^Jclca-, where ca is the 
concentration of actin monomers of size a [23]. This is valid when the persistence length of 
the chains is longer than the mesh size £. For a densely crosslinked gel, £ is also the typical 
distance between crosslinks, and therefore entanglement points: L e ~ £. In this case, 

£ max ~ ^ ~ — (ac A )- 1/2 (5) 

K K 

and 

G' ~ ^V 5 ~ ^{ac A f' 2 . (6) 
kT^ kT y ' v ; 

The precise dependence of the entanglement length on concentration in a solution of 
semiflexible chains is more subtle. We expect that L e may become substantially larger than 
£ for £ ^ £ p , since the transverse fluctuations of a semiflexible chain are greatly reduced over 
distances comparable to or smaller than the persistence length of the chain. For example, on 
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the scale of the mesh size £ < £ p , chains cannot form loops and knots [24-26]. We can estimate 
the entanglement length in the following way [22,24]. From the above energy in Eq. ([!]), the 
transverse fluctuations at temperature T of a chain confined (entangled) at one end grow 
as (L\) ~ kTL 3 /n, where L is the distance from the entanglement. Thus, the fluctuating 
chain segments of length L e between entanglements occupy a volume L e (L^) ~ kTL^/n. 
For a given concentration ca, the probability of an intersection with another chain is of order 
unity for L e ~ (k/UT) 1 ^ (ac^) -2 ^ 5 , which becomes larger than £ for £ <C £ p . Thus 

6 max ~ (kT/ K f 5 (ac A )- 2/5 (7) 

and 

G' ~ K(K/kT) 2/5 (a CA ) u/5 . (8) 

This model provides a consistent framework with which to understand the macroscopic 
viscoelasticity of biopolymer gels and solutions. Based on the semiflexible nature of several 
biopolymers, including F-actin, the model can explain both the large storage moduli as well 
as the observed strain hardening of networks at moderate to low strains [27], a feature in 
contrast with the behavior of flexible polymer networks [28]. For instance, at equal volume 
fractions, vimentin filaments, which are approximately an order of magnitude less stiff than 
F-actin, form networks with smaller shear moduli than F-actin, although vimentin networks 
can withstand approximately 10 times larger strains than F-actin networks before rupturing. 
Experimental observations of shear moduli and yield strain for varying actin concentration, 
as well as for modest changes in F-actin stiffness induced by binding of different nucleotides, 
are also in support of this model. 

This model makes several additional predictions that can be tested experimentally. First, 
as indicated above, for densely crosslinked gels, G' ~ tz 2 . Since it is now possible to measure 
k directly for actin and some other biopolymers by video microscopy [17], and there are 
a number of actin binding proteins and metabolites that can alter filament stiffness under 
conditions where filament length is held constant, the viscoelastic parameters can be directly 
measured as a function of k. Furthermore, the scaling behavior of entangled solutions and 
crosslinks gels as a function of concentration are predicted to differ. A third prediction is 
that the viscoelasticity of relatively dilute filament networks will be extremely sensitive to 
filament length even if the average filament length is much greater than the mesh size, and 
that this dependence will be greatest for the stiffest polymers. This is because for semi- 
flexible filaments the entanglement length required for effects on elasticity can be much 
greater than the mesh size, and this difference depends on k. Therefore, subtle changes 
in filament length can have large effects on viscoelasticity even when all filaments exhibit 
significant overlap. This feature may be one of the reasons the cytoskeletal actin filaments 
in cells are under the tight control of proteins that regulate their length. 
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FIGURES 

FIG. 1. Entangled network of semiflexible actin filaments. (A) In physiological conditions, 
individual monomeric actin proteins (G-actin) polymerize to form double-stranded helical filaments 
known as F-actin. These filaments exhibit a polydisperse length distribution of up to 70 /im in 
length. The persistence length of these filaments is of order 2 fj,m. (B) A dense solution (1.0 mg/ml) 
of actin filaments, approximately 0.03% of which have been labeled with rhodamine-phalloidin in 
order to visualize them by florescence microscopy. The average distance £ between chains in this 
figure is approximately 0.3 ^m. Note the nearly straight conformation of the filaments on this 
scale. 

FIG. 2. The plateau modulus G' of actin networks as a function of concentration in mg/ml [5]. 
The predicted scaling for entangled networks, from Eq. (@), is shown. In this case, 
nematic phase of actin filaments has been shown to form above a concentration of approximately 2 
mg/ml [19]. Our model is valid for the entangled isotropic regime [24], (kT/n) 2 /a < ca < kTj (no 2 ). 
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